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The emergence of new techniques for the fabrication of nematic droplets with nontrivial topology 
provides new routes for the assembly of responsive devices. Here we perform a numerical study 
of spherical nematic droplets on fibres. We analyse the equilibrium textures and find that, under 
certain conditions, the nematic can avoid the nucleation of topological defects. We consider in detail 
a homeotropic nematic droplet wrapped around a fibre with planar anchoring. We investigate the 
effect of an electric field on the texture of this droplet such type of system. In the presence of 
a DC field, the system undergoes a Freederickzs-like transition above a given threshold Ec. We 
also consider AC fields, at high and low frequencies, and find that the textures are similar to those 
observed for static fields, in contrast with recently reported experiments. 

PACS numbers: ... 


I. INTRODUCTION 

In the last decade, the structure and dynamics of ne¬ 
matic droplets have been studied, theoretically and ex¬ 
perimentally, due to their key role in novel applications 
based on polymer dispersed liquid crystals (PDLCs) [1]. 
The interplay between the spherical confinement and the 
molecular alignment at the surface of the droplets leads 
to the nucleation of topological defects. In general, two 
types of surface alignment (anchoring) are considered: 
homeotropic (perpendicular) and planar (parallel). Ne¬ 
matic liquid crystals confined to homeotropic spherical 
drops can exhibit a radial configuration with a point¬ 
like defect, of topological charge +1, at the center of 
the drop. The radial configuration is observed in large 
droplets, with radii R > 3fim [2]. For smaller droplets, 
however, the nematic texture adopts what is known as 
an axial configuration with a topological defect ring of 
cross-sectional winding number +1/2. The competition 
between the splay energy and the anchoring energy makes 
the point defect energetically unfavourable against the 
ring for sufficiently small drops miij. This transition 
may also be driven by applying an electric field [2] result¬ 
ing in distinct optical responses as the field is switched 
on or off. 

On the other hand, if the surface of the nematic drop 
induces parallel alignment, the droplet exhibits a bipo¬ 
lar configuration with two surface defects (boojurns) nu¬ 
cleated at antipodal points. In ellipsoidal droplets the 
defects are forced to the regions of higher curvature to 
relieve splay and bend elastic stress [5]. For some liq¬ 
uid crystals, such as 8CB, planar nematic droplets with 
a bipolar configuration can undergo a twist transition 
where the nematic director twists from one surface de¬ 
fect to the other. This transition is controlled by tem¬ 
perature, which turns twist deformations energetically 
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favourable against bend and splay deformations [6]. The 
twist transition does not affect the type of defects nucle¬ 
ated at the droplet surface. However, in systems where 
it is energetically preferable to have bend deformations 
rather than splay deformations, the nematic takes a con¬ 
centric configuration with a defect line connecting the 
poles of the droplet along the symmetry axis [7]. Pla¬ 
nar nematic droplets with bipolar configurations are the 
most common. Of particular interest are those with pro¬ 
late shapes that allow the switching between two distinct 
nematic configurations with the application and removal 
of an external electric field [1]. 

With the development of new soft-lithographic tech¬ 
niques it is now possible to produce in a controlled fashion 
new PDLC matrices with nematic droplets of prescribed 
shapes and non-trivial topologies. This gave origin to 
what was coined as topological PDLCs (TPDLCs) [8]. 
Of particular interest are toroidal shapes, for which the 
Hopf-Poincare theorem dictates that the net charge of 
the defects is Qi = 0, for homogeneous surface align¬ 
ment. It turns out that for a torus on the micrometer 
scale defects may be present in pairs of opposite charge 
[8]. However, as the system size increases the presence of 
topological defects becomes energetically unstable, and 
for a torus on the millimeter scale, the nematic texture 
is completely free from topological defects [9] . 

A topology which is also toroidal is that of a droplet 
adsorbed around a microfibre mm- In this case, how¬ 
ever, the system has two confining surfaces that can dif¬ 
fer in anchoring: (i) that of the droplet itself and (ii) the 
surface of the microfibre. Fabrication of these systems 
is easily achieved but control of the size of the droplets 
is limited. They were reported in m, where 5CB liq¬ 
uid crystal droplets were deposited on 1 fim thick fibres 
through aerosol evaporation. The droplets, with diame¬ 
ters between 10 to 30 jam, are in contact with air and the 
nematic texture adopts the axial configuration, with a 
ring defect around the fibre, indicating that the two sur¬ 
faces have antagonistic anchorings. At the nematic-air 
interface the 5CB molecules are aligned perpendicular to 
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the interface, implying that the surface of the fibre im¬ 
poses planar (or planar degenerate) alignment. The ad¬ 
vantage of this type of systems arises from the possibility 
of changing the surface alignment through the adsorption 
of certain molecules. For example, endotoxins are known 
to trigger ordering transitions in nematic droplets from 
the bipolar (planar alignment at the surface) to the radial 
(perpendicular alignment) configuration [ 12 ]. Also, the 
defect configuration may be strongly deformed with ap¬ 
plied electric fields due to the presence of the fibre, which 
pierces the droplet. The application of electric fields may 
change the configuration of the defect in isolated droplets 
but in most cases a rotation of the nematic texture has 
been reported m- 

Here we investigate the equilibrium textures of nematic 
droplets adsorbed at fibres for different boundary condi¬ 
tions. We show that the nucleation of topological defects 
in the LC matrix depends on the matching (or mismatch¬ 
ing) anchoring conditions at the fibre and surface of the 
droplet, and it may be avoided by allowing the nematic 
molecules to twist around the fibre. We will focus in par¬ 
ticular on the case of homeotropic droplets pierced by a 
planar fibre and investigate how the system responds to 
applied external electric fields. 


II. MODEL 

The nematic order of the fluid is expressed using 
a traceless, symmetric, tensorial order parameter [T4| 
called the Q-tensor, 

Qa/S 2^ ^OLf3^ T 2^ 5 (1) 

where n is a unit vector denoting the director, and S is 
the degree of nematic ordering. In some circumstances 
there may be biaxial ordering of degree H, with m and 1 
forming an ortho normal set with n. 

The free energy of the system is given by the functional 
over the fluid region IZ and the boundaries W 

:f = A 0a 
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where A, Li and L 2 , and a are positive coefficients for 
bulk, elastic and anchoring free energies respectively, r* 
is a reduced temperature, such that the nematic phase 
with S = S'nem IS favourcd when r* < 1 , and the isotropic 
(unordered) phase with S' = 0 is favoured when r* > 1 . 

is the value of the Q-tensor preferred by the anchor¬ 
ing at the boundary. In all cases reported here we choose 
the elastic constants such that L 2 = 2Li. This corre¬ 
sponds to liquid crystal compounds with a twist elastic 


constant lower than splay and bend elastic constants, as 
5CB. 

Within 7^, the fluid density p, velocity and Q-tensor 
evolve over time t according to the continuity, Navier- 
Stokes, and Beris-Edwards m equations: 


dtP + d0{pU0) =0 (3) 
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where p = p/3 is the isotropic fluid pressure, p is the 
dynamic viscosity, and T is the mobility of the nematic 
order. C is a dynamical parameter, dependent on the 
molecular details of the liquid crystal, which determines 
how the nematic orientation couples to shear. If C < 1, 
then the director will tumble indefinitely in the shear, 
while if C > 1 (which we shall consider in this paper) 
permits a bulk, state-state orientation of the director rel¬ 
ative to the shear m- 

At the boundaries W, non-slip and anchoring condi¬ 
tions apply. 


= 0, (7) 

iT^'y^'yQaI3 T L/2l^a^'yQI3j ^ t (S) 

where v is the inward normal to the substrate. As an 
alternative to Eqn. we may impose the boundary con¬ 
dition Q = This is equivalent to setting a = oo. 


III. SIMULATION METHOD 

We simulate the dynamics of the fluid by discretising 
space and time - the former into a cubic grid of nodes - 
and maintain p, u and Q as continuous quantities. We 
utilise a hybrid method in which Eqns. and are iter¬ 
ated using a lattice Boltzmann method and Eqn. [fusing 
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a a finite-difference method m, a method that has been 
previously used by our group when considering a nematic 
liquid crystal in contact with a substrate patterned with 
rectangular grooves that may fill without the occurrence 
of complete wetting m- In that study our numerical 
code was used to analyse the dynamical response of the 
system to an externally-applied electric field so as to iden¬ 
tify switching transitions between these filled states and 
was also applied by our group to study the effect of an¬ 
choring strength in the flow of a nematic liquid crystal 

m- 

To incorporate an electric field E into the simulation 
we take into account the electric contribution to the free 
energy 

-^elec ~ T ^A^nemQ^ (9) 

where eo is the electric permittivity of vacuum, e/ the 
relative permittivity and Ae the dielectric anisotropy. 
Ae > 0 corresponds to systems where particles tend to 
orient along the direction of the electric held, whereas for 
Ae < 0 particles orient perpendicular to the held. In this 
study, which is meant to model 5CB, we consider that 
particles align in the direction of the electric held with 
Ae = 11.5. An electric held E (AC) is defined as 

E{t) = £;mag cos(27r/t) (10) 

with Emag the magnitude of the electric held and / its 
frequency. Setting / = 0 reverts to a DC held of mag¬ 
nitude Emag- We follow experimental analyses and also 
many other numerical studies by assuming uniform E 
throughout the simulation box. 

Simulation parameters must be mapped to physical 
quantities of 5CB by an appropriate choice of scaling. Us¬ 
ing reported data for 5CB, A ^ 10^Nm“^, Li ^ 10“^^N 
with 1/2 = 2Li 10 ^m, nematic-air surface ten¬ 

sion 7 = 0.038Nm“^ and isotropic-nematic transition 
between 18 and 34° Celsius, we set our simulation pa¬ 
rameters to A = 0.188, Li = 0.085, L 2 = 2I/i, T = 0.25, 
(a = 1, p = 80, (with deviations due to the effects of com¬ 
pressibility, which are small), C = 1.5 (aligning regime) 
and 5'nem = 1- Hence the physical length-scale corre¬ 
sponding to a lattice spacing, is Ax ^ 10“^ m, and the 
time scale is At ^ 10“^ s. E^ag and / are parameters of 
interest that we shall vary. 

IV. NEMATIC DROPLETS ADSORBED AT 
FIBRES 

A. In the absence of E fields 

The texture or equilibrium configuration of a nematic 
droplet pierced by a fibre depends on (i) the anchoring 
conditions at both the fibre and the nematic-air inter¬ 
face, and (ii) the ratio between the size of the droplet 





FIG. 1: (color online) Nematic droplets at fibres with different 
anchoring conditions, a) Homeotropic alignment at both the 
fibre and the nematic-air interface; b) homeotropic alignment 
at the fibre and planar degenerate at the nematic-air interface; 
c) and d) planar degenerate at both the fibre and the nematic- 
air interface. The bars represent the director field n and the 
color code indicates the alignment of n with the orientation 
perpendicular to the fibre; red if n is perpendicular and blue if 
n is parallel to the fibre. The topological defect is represented 
by the isosurface in green. 

and the diameter of the fibre that pierces it, ddropidfibre- 
We start by considering different types of anchoring con¬ 
ditions. As discussed previously, nematic droplets with 
either homeotropic (perpendicular) or planar (parallel) 
alignment at their surface have topological charge of +1 
which induces the nucleation of defects either in bulk or 
at the surface, respectively. By piercing the droplet with 
a fibre, the structure of the defects can change or, in some 
cases, they can disappear. This depends on the type of 
anchoring at the fibre compared to the anchoring at the 
nematic-air interface. 

In Figwe show three nematic droplets adsorbed at fi¬ 
bres for different anchoring conditions. When the surface 
of the droplet induces homeotropic anchoring a defect is 
nucleated in bulk at the center of the drop. As discussed 
before, this can be either a point defect, which typically 
appears for large droplets and the nematic has a radial 
configuration, or a ring defect with a nematic texture 
that is axial. In Fig. [1]2 a fibre with homeotropic anchor¬ 
ing goes through the center of the nematic drop. The 
mere presence of a fibre going through the centre of the 
drop makes it physically impossible for the defect to nu¬ 
cleate in that same region. In the case of a homeotropic 
fibre this is compensated by the nucleation of two sur¬ 
face defects at the contact lines (see Figj^in green), and 
no defect is nucleated in bulk. Although we did not ex¬ 
plore it further, as our droplet is assumed to be rigid, 
in this situation its shape is expected to elongate. Such 
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FIG. 2: (color online) Homeotropic nematic droplet at a pla¬ 
nar fibre, a) Equilibrium configuration with a ring defect 
surrounding the fibre at the middle of the droplet. If the an¬ 
choring on the fibre is planar degenerate, b) there is no topo¬ 
logical defect in the nematic and c) the liquid crystal texture 
assumes an escaped configuration close to the fibre. The bars 
represent the director field n and the color code indicates the 
alignment of n with the orientation perpendicular to the fi¬ 
bre; red if n is perpendicular and blue if n is parallel to the 
fibre. The topological defect is represented by the isosurface 
in green. 

elongation would be balanced by the surface tension. 

Nematic droplets with planar degenerate anchoring 
can either exhibit a bipolar or concentric texture, as men¬ 
tioned before. However, in the presence of a homeotropic 
fibre, the defects of either configurations are no longer 
nucleated. In Figj^ we see that instead of nucleating 
topological defects, the nematic assumes an escaped con¬ 
figuration around the fibre, where the director field seems 
to rotate around the fibre. This is possible due to the de¬ 
generate orientation at the surface of the droplets. A 
similar effect is seen when one considers a fibre with pla¬ 
nar (degenerate) anchoring, shown in Figsj^ and d. In 
this case the director field rotates close to the contact 
lines, clearly seen in Fig{^, where the director assumes 
a perpendicular orientation to the fibre (symmetry) axis, 
while planar to its surface, and as it moves away from 
the contact line towards the center, the director aligns 
parallel to the symmetry axis. 

Finally, if a fibre with planar anchoring pierces a 
homeotropic droplet the nematic nucleates a topological 
ring defect close to the fibre as seen in Fig|^. This re¬ 
sembles the axial configuration of nematic drops, where 



FIG. 3: (color online) Free energy T of a homeotropic ne¬ 
matic droplet adsorbed at a fibre as a function of the applied 
DG electric field E. Fq — F[E — 0) is the free energy in the 
absence of the external field and Vti is the LG drop volume. 
Left: fixed system size ddrop = 0.4 /xm and different anchoring 
strengths o = 1 (□, red), 10 (o, black). Right: fixed anchor¬ 
ing strength a = 1 and different system sizes ddrop = 0.4/xm 
(dash line, red), 0.8/xm (dash-dot line, blue). The inset il¬ 
lustrates the dependence of Ec, given by the intersection of 
the same-coloured curves, on system size, as explained in the 
text. 


the fibre plays the role of the axially uniform nematic. 
This configuration is stable for any size of the droplet if 
the anchoring on the fibre is non-degenerate. However, 
if the fibre has planar degenerate anchoring the ring de¬ 
fect is unstable for large systems. It turns out that by 
increasing the effective distance between the surface of 
the nematic drop and the surface of the fibre, while fix¬ 
ing ddrop/dfibres the ucmatic molecules are able to ro¬ 
tate around the fibre at the center of the drop, thus re¬ 
moving the topological defect (Figj^) and assuming an 
escaped configuration that induces bend-splay deforma¬ 
tions (FigJ^). Such configuration is indeed not reported 
in the experiments m which hints that the fibre may in¬ 
duce a preferred planar orientation rather than degener¬ 
ate. For small systems with planar degenerate anchoring 
on the fibre the escaped configuration is metastable and 
can be achieved through the application of an external 
electric field. We note that if the surface of the nematic 
drop is contaminated by an agent that induces planar 
anchoring the liquid crystal texture would switch from a 
configuration with a ring defect to a configuration with¬ 
out any topological defects, inducing a different optical 
response. 


B. In the presence of DC E fields 

Homeotropic nematic droplets on fibres with planar 
anchoring are feasible from the experimental point of 
view. The appearance of a ring defect around the fi- 
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FIG. 4: (color online) Cross section of the liquid crystal texture of a homeotropic nematic droplet adsorbed at a fibre a) at 
equilibrium with E = 0, b) and c) under the influence of a DC electric held for held strengths respectively below and above the 
Freederickzs transition, and d) to /) evolving towards the equilibrium conhguration after switching off the electric held. The 
snapshots for the off sequence were taken for d) t = 5 ms, e) t = 10 ms, and /) t = 25 ms. The bars represent the director held 
n and the color code indicates the alignment of n with the applied electric held; red if n || E and blue if n T E. The topological 
defect is represented by the isosurface in green. 


bre enhances the potential of such systems for applica¬ 
tions, as an external held can induce an orientational 
state with an optical response that is signihcantly differ¬ 
ent from that of the ground state. Here we are interested 
on the response to applied electric helds. For comparison 
we start by assessing the effect of electric helds produced 
by direct current (DC) electrostatic potentials. Alter¬ 
nating current (AC) helds will be addressed in the next 
section. 

Here we consider that the electric held E is perpen¬ 
dicular to the hbre and couples with the nematic LC via 
the free energy density given in Eqj^ When applying 
an electric held the system will respond in accordance 
with the competition between the elastic energy, that 
acts as a restoring force, and the held energy. If the held 
is strong enough it will deform the orientational ground 
state and the system will take a characteristic time to 
reach the new state that is a product of such compet¬ 
ing mechanisms. Figure shows how the free energy F 
changes as the applied DC electric held E is increased. 
As expected, the system undergoes a Freederickzs type 
transition, well known for hat liquid crystal cells [TH [20] . 
Here, although the geometry of the system is more com¬ 


plex, the response to electric helds seems to be very sim¬ 
ilar. For weak held strengths, little reorientation of the 
nematic director occurs, the elastic strength is strong and 
the systems remains very close to its ground state conhg¬ 
uration. Above a critical held strength, Ec, the electric 
held aligning forces are stronger than the elastic ones 
and the nematic tends to align with the E held, thus 
deforming completely the zero held state. In Fig. [^we 
show the free energy F as a function of the applied elec¬ 
tric held E for (left) hxed system size ddrop = 0.4/im 
and different anchoring strengths a = 1, 10, and (right) 
hxed anchoring strength a = 1 and different system sizes 
ddrop = 0.4, 0.8/im. Fq is the free energy for F = 0 and 
Vjz is the volume of the nematic liquid crystal. For the 
same system size, increasing the anchoring strength re¬ 
sults in an increase of the transition held. For example, 
for ddrop = 0.4 /im increasing the anchoring from a = 1 
to a = 10 induces the transition held to change from 
Ec = 0.40 to Ec = 0.45. On the other hand, for the 
same anchoring strength a = 1, changing the system size 
from ddrop = 0.4/im to ddrop = 0.8/im has the reverse 
effect, and the transition held decreases from Ec = 0.40 
to Ec = 0.34. As the anchoring strength is increased 
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FIG. 5: (color online) Cross section of the liquid crystal texture of a homeotropic nematic droplet adsorbed at a fibre under the 
influence of an AC electric field for field strengths a) below the Freederickzs transition, b) close to and above the Freederickzs 
transition, and c) far and above the Freederickzs transition, and for frequencies z) / = 0 Hz (DC), ii) f ~ 100 Hz, in) / = 1 kHz, 
and iv) / = 100 kHz. The bars represent the director field n and the color code indicates the alignment of n with the applied 
electric held; red if n || E and blue if n T E. The topological defect is represented by the isosurface in green. 


the aligning effect of the surfaces gets stronger and it 
increases the energy cost of uniformly aligning the liq¬ 
uid crystal molecules along the held, which results in a 
higher transition held. Enlarging the separation between 
the conhning surfaces weakens this effect and the transi¬ 
tion held is lower. 

To address the response of the system to the electric 
held, we show in Fig. what happens to the nematic 
conhguration when the held is turned up to a high held 
strength and then switched off. When there is no elec¬ 
tric held and the system is at equilibrium, the nematic 
exhibits the axial conhguration of Figj^ with a defect 
ring around the hbre. As the electric held is switched 
on up to a value that is below the Freederickzs transi¬ 
tion {E < Ec) the LC molecules start to align with the 
held. This process is counteracted by the elastic stress 
and the system remains in the axial conhguration. How¬ 
ever, we do observe a contraction of the ring defect as an 
effect of the electric held (FigJ^). If the held strength is 
higher than Ec the nematic molecules mostly align with 
the held. In a nematic droplet, as the contraction of the 
defect becomes unstable, this would lead to a rotation 


of the defect ring, which would also expand towards the 
surface of the drop. Here, given the presence of the h- 
bre, the ring rotates (twists) in two different directions 
while expanding towards the surface of the drop, trans¬ 
forming the ring into a hgure-of-eight, as represented in 
Figj^. When switching off the electric held the elastic 
stress acts in order to restore the system to a state with 
lower elastic deformations. This results in the untwisting 
of the hgure-of-eight defect (FigJ^) and the contraction 
towards the hbre (FigJ^ and /). 

The time it takes for the system to achieve the new 
state depends on the processes involved. Always present 
is the elastic stress that acts as a restoring force that 
takes the system into a state of lower elastic deforma¬ 
tion. When applying an electric held the response of 
the system depends not just on the restoring effect of 
the elastic stresses but also on the perturbing/aligning 
effect of the electric held, and has a characteristic time 
given by Ton = 1/ [El’ where and te are 

respectively the characteristic times associated with the 
elastic and electric stresses. For very large held strengths 
^on ^ When the electric held is switched off the re- 
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laxation time depends solely on the characteristic time 
due to the elastic stress, Toff ^ tat- 

The relevant characteristic times can be estimated by 
the dimensional analysis of Eq. This indicates that 
the characteristic time associated with the elastic stress 
scales with tn oc A^/ (TLi), where A is a characteristic 
length. For (transient) configurations close to the (equi¬ 
librium) axial configuration A ^ dfibre- In the presence 
of an electric field E > the configuration is extremely 
different from the equilibrium one given that the defect 
expands to the surface of the drop. In spherical, or el¬ 
lipsoidal, nematic drops the defect would expand main¬ 
taining its ring shape and in that case A ^ ddrop- In 
the presence of a piercing fibre this is no longer the case 
and A is not as trivial to estimate. There are, how¬ 
ever, at least two characteristic time scales associated 
with elastic stresses. One that characterizes the time 
it takes for the distorted ring to regain its equilibrium 
shape, tati cx: (TLi), and another that is the time 

it takes for the figure-of-eight defect to untwist, tats- We 
have performed measurements of tat^ and r]sf 2 5 and found 
that, typically, = S.Otati and that for larger systems 
{ddrop ^ 1 /im) tati = 2.0 ms. We can also estimate the 
response time of the nematic to the applied E field to 
scale with te oc (TeoAeE^) . Our measured te con¬ 
firms the E~^ dependence. As E increases te becomes 
smaller and the effect of the electric field dominates over 
the LC elasticity. 


C. The effect of AC fields 

In LC displays it is possible to use either DC or AC 
fields. However, AC fields have a practical advantage 
over DC ones. AC fields typically produce more stable 
patterns, lower the threshold voltages and minimize elec¬ 
trochemical deterioration [21]. In the case of nematic 
droplets adsorbed at fibres, and in contact with air, the 
application of a DC field can disrupt completely the drop, 
which is pulled towards one of the electrodes. This is due 
to ionic impurities in the LC compound that aid the for¬ 
mation of space charges. However, the effect of free ions 
can be significantly minimized by the use of AC fields. 
Here we will not explore the effect of AC r’s DC fields on 
free ions and how these influence the LC texture, which 
will be done elsewhere. Instead, in this section we aim to 
address how AC electric fields affect the nematic textures 
presented in the previous section. 

We again consider that the electric field is applied per¬ 
pendicular to the fibre. We note that, for both DC and 
AC fields, if the electric field is oblique to the fibre, this 
results in a rotation of the defect. For simplicity we as¬ 
sume that the electric field is sinusoidal and is given by 
Eq{^(with no DC component), such that for frequency 
/ = 0 we regain the DC limit. Our study reveals that 
there are three types field strengths Emag to consider: 
i) below the Freederickzs transition field Ej^ag < E) 
just above but close to Ec^ and in) Emag > Ec. Ex¬ 


amples of these cases are show in Fig. for different 
field frequencies / = 0, 100, 1000,10000 Hz. For field 
strengths below Ec (FigJ^ ^), the nematic aligns weakly 
with the applied field E. This can result in a small de¬ 
formation/contraction of the defect ring. In this situa¬ 
tion the response time te and the restoring time tat^ are 
comparable when E = Emag- This means that as the 
frequency of the field is increased (FigJ^i^ toj^^-f;) the 
elastic stress dominates and the deformation/contraction 
of the defect is weaker. If, however, the field strength is 
strong {Emag > Ec) the nematic strongly aligns with the 
electric field and the defect expands (and twists) to the 
surface of the drop. In this situation the response time 
Te is smaller than the restoring time r]sf 2 - As the fre¬ 
quency is increased the electric field will spend a certain 
amount of time below Ec which results in an effective 
response time that is comparable to the restoring time 
r]sf 2 - As a consequence the defect shrinks and assumes a 
shape normally observed for DC fields close to the Freed¬ 
erickzs transition. The most striking phenomena occurs 
for fields above the transition but comparable to Ec. As 
the frequency is increased the defect regains its nearly 
circular shape for Emag < indicating that restoring 
time tati is effectively smaller than the response time te^ 
the elastic stress dominates, and the Freederickzs transi¬ 
tions is shifted. 


V. CONCLUSIONS 

Nematic droplets suspended on fibres are new LC sys¬ 
tems that very much resemble water droplets on cob¬ 
webs. They differ from the nematic droplets dispersed 
in a polymer matrix, as the fibre piercing through the 
droplet confers the droplets a torus-like topology. Here 
we have studied the equilibrium textures of spherical ne¬ 
matic droplets on fibres under several anchoring condi¬ 
tions. We have shown that if a nematic droplet with 
homeotropic anchoring is pierced by a homeotropic fi¬ 
bre two surface defects are nucleated on antipodal posi¬ 
tions and at the contact lines, which in principle could 
be removed by elongating the particle. However, if the 
droplet has planar degenerate anchoring, and the fibre 
has either homeotropic or planar degenerate anchorings, 
the nematic can avoid the nucleation of topological de¬ 
fects by twisting along the symmetry axis. Finally, if a 
homeotropic nematic droplet is pierced by a planar fibre 
a ring defect is nucleated close to the fibre in the centre of 
the drop. Such a defect may be avoided if the fibre has 
degenerate planar anchoring. In this case the nematic 
will rotate around the fibre and the orientational field 
assumes a smooth escaped configuration. 

We have also studied the effect of applying an electric 
field to a homeotropic nematic droplet pierced by a pla¬ 
nar fibre. We have started by considering the application 
of a DC field, perpendicular to the fibre, and have seen 
that the system undergoes a Freederickzs like transition 
at a field strength Ec that depends on the size of the sys- 


tern and on the anchoring strength at both the surface of 
the droplet and of the fibre. Below Ec the defect ring un¬ 
dergoes a small deformation/contraction. Above Ec the 
ring expands towards the surface of the droplet, similar 
to what is observer for spherical nematic homeotropic 
droplets. However, the presence of the fibre forces the 
ring to twist and assume a figure-of-eight shape. 

Finally, we have considered the effect of an AC field. 
In our model, the frequency of the AC field seems to 
play a small role. Its effects are more pronounced at field 
strengths just above the Freederickzs transition, where 
instead of the figure-of-eight defect spanning throughout 
the system we observe a ring defect (untwisted) indicat¬ 
ing that the AC field shifts the transition to higher field 
strengths. 

We note that in our simulations we did not observe 
dynamical phenomena similar to that reported in HO]. 
Ceng and co-workers applied to the system, modelled in 
this work, an AC field of fixed frequency / = 50 Hz. They 
reported that for field strengths below a critical value the 
ring defect shifted along the fibre in one direction and 
on the other a spacial soliton seemed to appear rotating 
with a frequency that depended on the field strength. 
We can suggest (at least) two reasons for the discrepancy 
between the dynamics reported for the real system and 
the behaviour of our model. 

The first is related to the size of the systems used in our 


simulations, which are two orders of magnitude smaller 
than the physical systems. The size of the simulated 
system is limited by the resolution required to describe 
the defect, of the order of the bulk correlation length, 
^ ^ 15 nm, which is orders of magnitude smaller than 
the droplet radius, in the range 10 to 30 /rm. 

The second reason is related to the presence of ionic 
impurities in the real systems. Ionic impurities may be 
present in LCs due to molecular degradation, charge in¬ 
jection, or industrial processes. They may be introduced 
into the LC matrix during sample assembly through the 
contamination of substrates or of polymer alignment lay¬ 
ers [22]. In the presence of electric fields free ions will 
induce flows [23l [24] and local distortions of the net elec¬ 
tric field, which may dominate the dynamics. Including 
the effect of ions in our simulations is a major task, well 
beyond the scope of this work. 
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